Optimal unambiguous discrimination of two subspaces as a case in mixed state 

discrimination 
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We show how to optimally unambiguously discriminate between two subspaces of a Hilbert space. 
In particular we suppose that we are given a quantum system in either the state where \ipi) 

can be any state in the subspace Si, or \ip2), where \ip2) can be any state in the subspace S2, and 
our task is to determine in which of the subspaces the state of our quantum system lies. We do 
not want to make any error, which means that our procedure will sometimes fail if the subspaces 
are not orthogonal. This is a special case of the unambiguous discrimination of mixed states. We 
present the POVM that solves this problem and several applications of this procedure, including 
the discrimination of multipartite states without classical communication. 

PACS numbers: 03.67.-a,03.65.Bz,42.50.-p 



I. INTRODUCTION 



The discrimination of quantum states is an area that 
has received considerable attention. For a recent re- 
view see pj, for example. Pure state discrimination has 
proven to be useful for both quantum cryptography and 
for quantum algorithms 0, . The basic protocol is the 
following. One has a list of possible quantum states, and 
one is presented with a system that is guaranteed to be 
in one of them. Our task is to determine which. If the 
states are not orthogonal, they cannot be discriminated 
perfectly, and it is necessary to specify what kind of strat- 
egy we wish to use. One possibility is minimum error 
discrimination p|. In this procedure, we always specify 
a state, but, because the states are not orthogonal, there 
is some chance that we will make a mistake. Minimum- 
error procedures minimize the chance of making a mis- 
take. Another alternative is unambiguous discrimination. 
This procedure never makes a mistake, but it can some- 
times fail. If the list of states contains N possibilities, 
then this procedure has N + 1 possible outputs, one for 
each of the states on list plus a failure output. If we re- 
ceive an output that specifies one of the states, then we 
know what the state of our system was, and if we receive 
the failure output, then we have no idea what it was. An 
optimum unambiguous procedure is one that minimizes 
the probability of receiving the failure output. We shall 
be considering unambiguous discrimination here. 

The problem of optimal unambiguous discrimination 
between two pure quantum states was solved by Ivanovic, 
Dieks and Peres 0, IE 0- They assumed that each of 
the two states were equally probable. The case in which 
they are not was treated by Jaeger and Shimonyj8[ (for a 
somewhat simpler derivation of their result see [9j)- The 
case of more than two pure states is not simple and has 
been considered by a number of authors. There are a few 
general results, but explicit procedures are available only 



for special cases [lJJ, ^| . One important result is that for 
unambiguous discrimination to be possible, the states on 
the list must be linearly independent |PJ. There are also 
lower bounds on the failure probability |l3| . 

The discrimination of two mixed states has only been 
considered more recently. Before we cite the previous re- 
sults it will be useful to introduce some terminology at 
this point. The support Sj of the density operator pj, 
describing a quantum state, is the subspace of the entire 
Hilbert space Tt spanned by the eigenvectors of pj belong- 
ing to nonzero eigenvalues (for j = 1,2). The rank of the 
density operator is equal to the dimension of its support. 
The subspace orthogonal to its support, Sj, is the kernel 
of the density operator pj such that TL = Sj Sj. We 
shall denote the projector onto Sj by Pj and the pro- 
jector onto Sj by Pj. Equipped with these definitions 
we can now make the following general statement. For 
two mixed quantum states unambiguous discrimination 
is possible with a finite probability of success if and only 
if the supports of their density operators are not identi- 
cal. Indeed, in such a case the kernel of at least one of 
them is not empty and a projective measurement along 
this kernel unambiguously identifies the other state. On 
the other hand, if the supports are identical, then so are 
their kernels, and there is no direction in TL that could 
unambiguously identify at least one of the density oper- 
ators. 

We can now return to a listing of earlier results. The 
POVM for unambiguously discriminating between a pure 
state and a rank two mixed state was derived in |14| and 
subsequently generalized to the case of unambiguously 
discriminating a pure state from any mixed state |3|]. 
Lower bounds on the failure probability for the unam- 
biguous discrimination of two mixed states were derived 
by Rudolph et al. [15J and for an arbitrary number of 
mixed states by Feng et al. ^(| • Raynal et al. proved 
two theorems that make it possible to reduce the problem 
of unambiguously discriminating between two arbitrary 
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mixed states of rank fci and k% to the discrimination of 
two states of the same rank, k < min^i, feg), in a 2k- 
dimensional space. Building on the results of 0] and 
[l6| , Herzog and Bergou |l8j] found explicit solutions for 
some special cases along with necessary conditions for the 
saturation of the lower bound. In particular, these results 
showed that unlike in the case of two pure states whether 
or not the lower bound can be attained depends not only 
on the value of the prior probabilty of the states but also 
on their structure. There are mixed states for which the 
lower bound can not be reached for any value of the prior 
probability. In 19] the optimal measurement operators 
were constructed explicitly for some special cases. 

In this paper we shall consider the unambiguous dis- 
crimination between two subspaces. What this means is 
the following. A state is chosen from one of two sub- 
spaces, and we wish to determine to which of the sub- 
spaces the state belongs. Within each subspace each 
state is equally likely, though one subspace may be more 
likely than the other. One place in which this type 
of problem has arisen is in the consideration of pro- 
grammable discriminators |20| . In this case, one is given 
three qubits, the first two are arbitrary but the third is 
guaranteed to be identical to either the first or the second 
qubit, and the problem is to determine which two qubits 
are identical. The problem can be solved by realizing that 
one is, in fact, discriminating between two subspaces, the 
first being the subspace of three-qubit states that is sym- 
metric in the first and third qubits, and the second being 
the subspace that is symmetric in the second and third 
qubits. 

Subspace discrimination is a special case of the dis- 
crimination of two mixed states; in this case the density 
matrices are just proportional to the projection operators 
onto the subspaces. Making use of the results of Raynal 
et al. ^3 j we can restrict our attention to the case of 
two subspaces of dimension k in a 2k dimensional space. 
In particular, let Si, and S2 be fc-dimensional subspaces 
of the entire Hilbert space, Ti, which has dimension 2k. 
We can assume that the intersection of Si and 62 is just 
the zero vector, a situation which we henceforth refer to 
as general position. We assume that pi = (l/k)Pi oc- 
curs with probability 77, and P2 — (l/fc)P2 occurs with 
probability 1 — 77, where Pj is the projection onto Sj, for 
j = 1,2. The POVM that distinguishes them has three 
elements, IIi, II2, and IIo = / — IT — II2, all of which 
are positive operators. The probability of identifying pj 
if we are given pj is pj — Ti(pjllj), and the probability 
of failing to identify it is qj — Tr(pjllo), for j = 1, 2. The 
condition that a state never be misidentified implies that 
H1P2 = II2/91 = 0. The average failure probability is 

Q = mi + (1 - r,)q 2 , (1.1) 

and our object is to find, for a given 77, a POVM that 
minimizes Q. From the results in and ^(|, we have 
that 

Q>2^r 1 (l-r 1 )F(pi,p 2 ) , (1.2) 



where the fidelity between the two density matrices is 
given by 

F(pi,p 2 )^Tr((p 1 1 /2 p 2 p 1 1 /2 ) 1 ^) . (1.3) 

The conditions, under which this bound can be satu- 
rated, have been investigated in 18] . In the opti- 
mum POVM has been given explicitly for certain spe- 
cial cases. As we shall see, unlike for the case of two 
pure states, this bound cannot al way s be reached for two 
mixed states, in agreement with [18J. The optimal mea- 
surement procedure depends on the value of 77, a feature 
it has in common with the procedure for discriminating 
two pure states. For two pure states there is always a 
range of 77 where equality in l|1.2|l can be reached. In 
the case of two mixed states, however, the optimal mea- 
surement procedure also depends on the structure of the 
two density operators, a distinctive feature that has no 
equivalent in the case of two pure states. For 77 near 
or 1 the optimal measurements are projective ones. In 
the intermediate regime, the optimal measurements are 
intermittently POVM's or projective measurements and, 
in general, their failure probability is higher than the fi- 
delity bound (|1.2fl . Only under very special conditions 
shall we find that the fidelity bound can be saturated. 

The main technical device that we shall use to find the 
optimal measurements is that of Jordan bases. These 
bases take the following form. The states |V>i), . . . \ipk) 
form an orthonormal basis for St, \ipk+i), ■ ■ ■ IV^fc) form 
an orthonormal basis for S 2 , and, in addition these states 
have the property that 

(ipi\ipk+j) = Sij cos8 t , (1.4) 

where cos6q > cos 8 2 > . . .cos 8k, and 1 < i, j < k. The 
states \ipt), \ip k ) and \ipk+i), \1p2k) are called Jor- 
dan bases and the angles 8i are called the Jordan angles. 
Bases satisfying these conditions can be constructed for 
any two subspaces j^- Note that the basis vectors \ipi) 
and \ipk+i) are eigenvectors of the operators P1P2P1 and 
P 2 PiP 2 , respectively, where 

WW = cos 2 8^) , 

P2PlP2^k +i ) = COS 2 ^|7/> fe+ . t > , (1.5) 

where 1 < i < k. 

The paper is organized as follows. Sec. [H]is devoted 
to the derivation of the general results. Besides giving 
the general theory for distinguishing two k dimensional 
subspaces in a 2k dimensional Hilbert space, these results 
also hold for a rather general class of density operators, 
so they are directly relevant to the problem of optimal 
unambiguous discrimination between two mixed states. 
In Sec. IIIII we present some possible applications of the 
results. Finally, in Sec. II VI we give a brief summary and 
outlook for future research. 



3 



II. DISTINGUISHING SUBSPACES AND ITS 
RELATION TO THE DISCRIMINATION OF 
MIXED STATES 

We shall actually solve a somewhat more general prob- 
lem than the discrimination of two subspaces. Let Si and 
S2 be k dimensional subspaces of a 2k dimensional com- 
plex Hilbert space Ti. which are in general position, as dis- 
cussed in the introduction, and let • • • , \tpk)} and 
{{tpk+i), ■ ■ * ) \4>2k)} be Jordan bases associated to these 
subspaces. Consider the two density matrices 

k 

pi = ^2ai\ipi)(ipi\ , 

i=l 
k 

P2 = 5^/3t|V<+fc)<&-+fc| , (2.1) 

i=l 

where a % > 0, J2i a * = 1 > A > °> and J2i A = 1- 
Clearly, p\ has support in Si and P2 has support in S%. In 
this case the orthonormal frames for the density matrices 
given by the spectral theorem coincide with the Jordan 
frames of the supports. If 

k 

pi = (i/*)£hfc>M 

i=l 
k 

P2 = (l/k)J2\^+k)(^+k\, (2.2) 

i=l 

we say that pi and pi are uniformly mixed states, and dis- 
criminating between uniformly mixed states corresponds 
to the case of discriminating between the subspaces S\ 
and S2. Inserting the weights ai and Pi will allow us 
to address several issues in the general theory of mixed 
state discrimination. 

We will now construct an optimal POVM to unam- 
biguously discriminate between pi and pi- The POVM 
elements are LTi, II2 and Uq = / — III — II2 with the 
properties as discussed in the introduction. Ill and II2 
are self adjoint operators supported on S2 and Si, re- 
spectively. In order for them to form an optimal POVM 
they must be positive and, crucially, the rank of ITo must 
not exceed k [17]]. We wish to choose II 1 and IT2 so that 
the expression P{q) — vfTv(Jlipi) + (1 — 77)Tr(n 2 p2) is 
maximized. This is the same as minimizing the average 
failure probability Q(rj) = r/Tr(IIoPi) + (1 - r/)Tr(n p 2 ) 
that was introduced in 

Let Ti be the linear subspace spanned by \ipi) and 
iV'i+fe). The Ti's, with 1 < i < k, are mutually orthogo- 
nal, two dimensional subspaces, which are invariant un- 
der both pi and pi- Let us define states \zi) and \yi) 
as 

\ipi) = sin 9i\zi) + cos 9i\ip i+k ) 
\ip i+ k) = sbx9i\yi)+ cos 9i\ipi), (2.3) 

where \zi) (\yi)) is the orthogonal complement of \ipi+k) 
in Ti. Furthermore ■ • • , |^fe)} forms an or- 



thonormal basis for S2, and {\yi), ■ ■ ■ , \yk}} forms an or- 
thonormal basis for Si . 

If we write ITi = J2i j=i a ij\ z i)( z j\ and ^-2 = 
£i.j=i hj\yi){yj\, then 

k 

pi = Ti(U lP i) = sin 2 9 l , (2.4) 

and 

k 

p 2 = Tr(n 2 p 2 ) = J2 sin2 6i ■ ( 2 - 5 ) 

These equations do not depend upon the off-diagonal 
terms of 111 and II2. If IT and II2 are to be elements 
of a POVM they must be positive so an, &« > is a min- 
imum requirement. The presence of off-diagonal elements 
imposes additional restrictions on the diagonal elements 
if we wish to ensure positivity. Since the off-diagonal ele- 
ments do not play a role in pi and pi , it suffices to search 
for our optimal POVM among the diagonal operators. 

Let pi = (i/> i |IIi|V'j) and p i+k = (^ i +fc|n 2 |V'i+fc) be 
the individual success probabilities of the Jordan basis 
states for 1 < i < k. Let qi = (i/^i | ITo I V'i } and q~i+ k = 
(ipi+klliol^i+k) be the individual failure probabilities of 
the Jordan basis states for 1 < i < k. Here we introduced 
the overbar notation in order to distinguish the partial 
success and failure probabilities, pi and q~i, of \ipi) from 
the total success and failure probabilities, pj and qj, of pj 
for i < k and j — 1,2. Obviously, pi + qi = 1 holds. We 
can then express the POVM operators as IT = Jj =1 II M 
and n 2 = £, fe =1 n 2jJ where 

U hi = (2.6) 
sin 0i 

and 

n 2 ,< = ^rlwKwl- (2.7) 

sm 0i 

We can also set n = Y^i=i^o,i, where n ,i = Jt 4 — 
Ili i — 112,2, and 1^ is the identity in Ti. 

We now need to determine the values of q~i and q~i+k- 
This can be done by noticing that what we have done is to 
reduce our problem to k problems of optimally discrimi- 
nating two vectors, in particular the vector \ipi) from the 
vector \i/ji+k) in Ti. In more detail the situation is the fol- 
lowing. In our overall ensemble, with pi occurring with 
probability 77 and p 2 with probability 1 — 77, the prob- 
ability of occurrence for \tpi) is rjoii and the probability 
for the occurrence of \tpi+k) is (1 — Therefore, the 

probability for the occurrence of a vector in Ti is just the 
sum of these probabilities, 

p(Ti) = r, ai + (1 - (2.8) 

Now, the probability that \ifti) occurs given that Ti has 
ocurred is p(i\Ti) = rjai/p(Ti) and the probability that 
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\i>i+k) occurs given that Tj has occurred is p(i + fc|Tj) = 
(1 — r/)Pi/p(Ti). Consequently, in T- t we want to un- 
ambiguously discriminate \tpi) occurring with probability 
p(i\Ti) and \ipi+k) occuring with probability p(i + k\Ti) 
so as to minimize the failure probability 



Qi(r]) = p(i\Ti)qi +p(i + k\Ti)q i+ k. 



(2.9) 



This problem was first solved by Jaeger and Shimony § , 
(a somewhat simpler solution is given in |9j), and we can 
now make use of that solution. Let 



Pi COS 2 0i 



Pi 



_a t + Pi cos 2 Qi Pi + a t cos 2 i 
and, in addition, let 



■qoti 



[*,%] (2.10) 



(2.11) 



For a given 77, the value of qi which minimizes Qi(r]) is 

[ 1 v < a 

Qiin) = < <lT(v) if V ^ in h (2-12) 



cos" 



77 > di 



and qk+i (v) — cos2 @i I Ii iv) ■ Furthermore the rank of 
LTq^ is one, with the nonzero eigenvalue given by 



A; 



cos 2 Qi 



2 cos" 64 



sin 2 0i 



(2.13) 



with the corresponding eigenstate 



cos6i(l-qi). 
\d) = t^t m+k) 



sin 



qi - cos-' Vj 
sin 2 6i 



\*k) ■ (2-14) 



This specifies the POVM within T t . 

The optimal overall failure probability can now be ex- 
pressed as 

Q o P t = j-QOPt p{Ti) 

i=l 
k 

= ^yiouqiirf) + (1 - v)fo9k+i(v)] , (2-15) 



where Q° pt is the failure probability that results when Eq. 
(|2.12|) is substituted into Eq. I|2.9|l . Its explicit expression 
is given by 



rjOi + (1 - rj)Pi cos 2 l if r\ < c 4 



Q° pt p{Ti) = { 2^(1-7?)^^! cos 0i\ if Cl < n < di • 
77a, cos 2 0i + (1 - 77)/% if 77 > 

(2-16) 

The center line is the geometric mean of the two terms in 
either the first or the last line and, therefore, represents 
an absolute minimum for Qi. We obtain the absolute 
possible minimum of the total failure probability if we 



sum the center lines for all 1 < i < k. The summation 
yields 



Q° pt = 2 v/t7(1-77) J2 V*Ji\ cos 0< 



(2.17) 



Clearly, this absolute minimum can only be realized if 
and only if the intersection of all of the intervals Ii is not 
empty and the operating value of 77 is in this intersecton. 

The interpretation of this expression is straightfor- 
ward. Making use of the structure of the Jordan bases, 
1 1 

we obtain p\pip1 = J2i a iPi cos 2 0i\ipi)(xpi\, after a sim- 
ple calculation. Comparing this expression with Eq. (|1 
tells us immediately that Eq. I|2.17|l can be cast to the 
form 



Q° pt = 2y/T,(l-T))F( Pl ,p2) 



(2.18) 



where F(p 1 ,p 2 ) = Y,i V^iAl cos0i|. Here F{p ll p 2 ) is 
the fidelity between p\ and p 2 , constructively proving 
that the fidelity bound, Eq. 1)1. 2f> . can be saturated. To 
obtain an explicit expression for the fidelity would be a 
hopeless task, in general. What made it possible here is 
the fact that the density operators are diagonal in the 
Jordan bases of their support and we could take full ad- 
vantage of the ensuing Jordan structure. Furthermore, 
the above expression for the optimum failure probability 
holds only if 77 is an element of the intersection of all of 
the Ii intervals, Iq — f] i=1 Ii- Io may be empty, it of- 
ten is. Note, however, that in the case where p\ and p 2 
are uniformly mixed states, which is the case of subspace 
discrimination, 



Ii = 



1 



COS 2 0i 



1 + COS 2 0i 1 + cos 2 I 



c I 



i+1 



(2.19) 



so Iq = I\ ^ and the fidelity result holds in the entire 
1 1 interval. 

For the case when k = 2 there are only two such inter- 
vals, 1 1 and I 2 , and we will now give a complete classifi- 
cation of their intersection pattern. To this end, we first 
introduce a one-parameter characterization of pi and p 2 - 
Let us set a% = a and, consequently, a 2 = 1 — a. Simi- 
larly, we set Pi = P and, consequently, p 2 = 1 — p. We 
can now introduce a two-dimensional parameter plane, 
ap, where the square in the first quadrant, bounded by 
< a < 1 and < P < 1, corresponds to physically ac- 
ceptable choices for mixed states. So our task is reduced 
to finding the regions within this square with qualita- 
tively different overlap patterns. When cos 2 0\ > cos 2 2 , 
which is the convention that we adopted at the begin- 
ning, the patterns can be sorted into five categories: 
i) d\ < c 2 , i. e. I\ is to the left of I 2 and their intersection 
is empty. This happens when 



- a cos 2 6>i cos 2 9 2 
P < Pi{a) = j- 5— - 

1 — a(l — cos 0i cos^ 9 2 ) 



(2.20) 



Pi (a) is a hyperbola in the aP plane and it is the divider 
between this region and the next, when 
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ii) ci < C2 < di < d%, i. e. the right end of I\ partially 
overlaps with the left end of I2 and their intersection is 
the overlap. This happens when 

f3<k(a) = — -r-. (2.21) 

cos 0i — a(cos^ 0i — cos 2 ) 

@2 (a) is a hyberbola in the a/3 plane and it is the divider 
between this region and the next, when 

iii) C2 < ci and d\ < d 2 , i- e. 1\ is inside I2 and the 
intersection coincides with 1\ . This happens when 

cos 2 + a(cos^ 0i — cos 2 ) 

/?3 (a) is a hyberbola in the a/3 plane and it is the divider 
between this region and the next, when 

iv) ci < C2 < g?2 < di, i. e. the left end of I\ partially 
overlaps with the right end of I2 and the intersection is 
the overlap. This happens when 



P < fa(a) 



a 

cos 2 0\ cos 2 62 + a(l — cos 2 0\ cos 2 



tf2 i- Ct(l — COS" (71 COS' C2J 

(2.23) 

/34(a) is a hyberbola in the a/3 plane and it is the divider 
between this region and the next, when, finally 
v) d 2 < c i 5 i- e. I\ is to the right of I2 and the intersection 
is empty. This happens when 



/3>/3 4 



(2.24) 



We note that when cos 2 0\ — cos 2 02 the two inner 
dividers, /3~2 and fe, both degenerate into the diagonal of 
the square, /3 = a. Our findings are summarized in Fig. 
^where the five regions of the parameter space, resulting 
from the four dividers fix,.. -,04, are displayed for the 
representative values cos 2 Ox = f and cos 2 6 2 = \- 

Next we give two illustrative examples for k = 2. In 
these examples we prescribe ITo by giving its eigenvectors 
and eigenvalues \ for the possible 77 interval config- 
urations. We adopt the notation \Q(r])) (Xi(i])) for the 
eigenvector (eigenvalue) corresponding to = cfi pt - 

In our first example we consider the case where px and 
P2 are uniformly mixed, corresponding to subspace dis- 
crimination. Then Ix = [ci,dx] C I 2 = [02,^2], where 



cos 2 0i 
1 + cos 2 0i 



1 + cos 2 0i 



(2.25) 



The interval < 77 < 1 is divided into five subintervals, 
[0, c 2 ], [c 2 ,ci], [ci,di], [di,d 2 ], and [^2,1]- Tablenjsum- 
marizes the behavior of IIo in each of the subintervals. 

For our second example we choose px and P2 so that 
dx < C2, which can be easily arranged by picking (3 < 
/3i(a). We then have intervals [0, ci], [ci,di], [^1,02], 
[c2,g?2]j and [c?2,l]- This is a situation in which the fi- 
delity bound for the failure probability, Eq. 11.2f) . can 
never be achieved. The behavior of il is summarized in 
Table Inl 
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FIG. 1: Regions of the parameter space. Dotted line: (3i(a), 
short dashed line: /32(a), medium dashed line: /33(a), long 
dashed line: /34(a). In the region below /3i(a) and above 
/34(a) the intersection of Ii and I2 is empty, the fidelity bound 
can not be reached; in the regions between /3i(a) and /32(a), 
and between /33(a) and /34(a) the intervals 7i and I2 partially 
overlap, the fidelity bound can be reached in these overlaps; 
and in the region between /32(a) and /33(a) the interval ii is 
inside I2, the fidelity bound can be reached in the entire I\. 
For the figure we used the values cos 61 = ^ and cos 82 = \ ■ 



TABLE I: Eigenvalues and eigenvectors of IIo in the various 
intervals of r\ for Example 1. Both states are uniformly mixed 
and the intersection of I\ and I2 is the entire I\. The fidelity 
bound can be reached in the intersection. 



V 


Ai 


A 2 


|Ci> 


l<2> 


[0,C2] 


1 


1 


iv>i> 


1^2) 


[C2,Cl] 


1 




IV>i> 


IC2(77)> 


[ci,di] 




Aa(ij) 


iCifo)) 


\Uv)) 


[di , d 2 ] 


1 


A 2 (r7) 




\C2(V)) 


[d 2 ,l] 


1 


1 


1^3) 


l^4> 



TABLE II: Eigenvalues and eigenvectors of IIo in the various 
intervals of r\ for the second example in the text. The inter- 
section of Ii and I2 is empty and the fidelity bound cannot 
be reached in this case. 



V 


Ai 


A 2 


Id) 


IC2> 


[0,cx] 


1 


1 




\>h) 


[cx,di] 


Mv) 


1 




1^2) 


[di,c 2 ] 


1 


1 


1^3) 


1^2) 


[C2 , d 2 ] 


1 


Mv) 


1^3) 




[da,l] 


1 


1 







The trend is clear from these two examples. For each of 
the five regions of the a/3 parameter plane the operating 
value of 77 can have five possibilities. It can be outside 



G 



of the intervals I\ and I2 (three such intervals if I\ and 
I 2 do not intersect and two if they do), it can be in the 
nonoverlapping regions of I\ and I2 (two such intervals) 
and, finally, it can be in the intersection of I\ and I2 
(zero or one such interval). The five parameter regions 
and the five possibilities for r\ in each of these regions give 
us altogether twenty five characteristically different cases. 
In only three of them can the fidelity bound be reached. 
In twelve cases the optimum measurement is a standard 
von Neumann projection and in the remaining ten cases 
it is a combination of projections in some dimensions and 
POVMs in the others. 

We believe that these trends are general and they hold 
for the discrimination of any Rank 2 mixed states not 
just the ones where the Jordan basis coincides with the 
spectral representation but it will be much harder to find 
explicitly the different regions in the parameter space and 
the different intervals of the prior probability ij. We also 
conjecture that for the discrimination of two Rank N 
mixed states there are (27V + 1) N cases altogether and 
the fidelity bound can be reached in only 2N — 1 of them. 
Since the growth in the number of possibilities as a func- 
tion of N is faster than exponential it seems as though 
it would be extremely difficult to give a complete classi- 
fication of the cases for N > 2. Furthermore, since the 
number of cases when the fidelity bound can be reached 
grows only linearly with N, the weight of the density 
operators for which the fidelity bound can be attained 
quickly becomes negligible with increasing N. 



III. APPLICATIONS 

Let us now consider a simple example that we will 
be able to use as the basis for applications of subspace 
discrimination. Let 7i be a four-dimensional space with 
the orthonormal basis {\j) \ j = 0, . . . , 3}. For the first 
subspace, Si, we choose the span of the vectors |0) and 
|1), and for the second, S% we choose the span of the 
vectors \u ) = (|0) + |2))/V2 and \ui) = (|1) + |3))/V2. 
The states {|0), |1)} and {|uo)i \ u if\ form Jordan bases 
for the subspaces Si and £2- Defining the vectors 



\Vi) = |2> life) = |3) 



| Zl ) = _(|0)-|2)) 



\z 2 ) = |3», (3.1) 



we have that Si is the span of \yi) and 1 2/2 } , and S2 is the 
span of \zi) and |z2). Application of the formulas in the 
previous section gives us the POVM for discriminating 
between S x and S 2 - The POVM exists for 1/3 < 77 < 2/3, 
where the detection operators are given by 



Hi = V2(V2-^^- 1 />, 



n 2 = V2^/2 



V 



1-r, 



Pi 



(3.2) 



The corresponding failure probability is Q — v/2r7(l — 77). 
We also mention here that this solution was already de- 
rived in 0| using a slightly less general approach. 

One possible application of this POVM is the follow- 
ing. Suppose that Alice and Bob cannot communicate 
with each other, but they can communicate with Char- 
lie. Charlie wants Alice and Bob to share a secure bit 
string. He sends to Alice and Bob one particle each from 
either of the two-particle states 

l*o} - -^(|0) a |l> 6 + |l>„|0) 6 ) 

|*i) = \{\u Q ) a \ui) b + \ui) a \u )b)- (3.3) 

If Alice and Bob both succeed in identifying which state 
was sent, they share a bit, \^o) corresponding to and 
l^i) corresponding to 1. The reduced density matrices 
that Alice and Bob must distinguish are po = (l/2)Pi, 
which results if l^o) is sent, and po = (1/2)P2, which 
results if l^i) is sent. The above POVM does this op- 
timally (we shall assume that 77 = 1/2). The procedure 
would be the following. Charlie sends one of the two 
states to Alice and Bob (one particle to each). They 
independently perform their measurements. They then 
tell Charlie whether they succeeded, and he tells each of 
them whether the bit is valid or not. The bit is valid 
when both Alice's and Bob's measurements succeeded, 
and invalid otherwise. 

The security comes from the fact that \^q) and are 
not orthogonal. An eavesdropper, Eve, cannot perfectly 
distinguish these two states. Her measurement proce- 
dure will either sometimes produce errors or sometimes 
fail. She must, however, send particles on to Alice and 
Bob. There is no state she can send them that will guar- 
antee that one or both of their measurements fail, so that 
sometimes Alice's and Bob's measurements will tell them 
that they have received a state different from the one that 
Charlie sent. By comparing some of their bits with those 
of Charlie, they can tell whether this has occurred. One 
possibility is that they can use some of the invalid bits, 
in particular the ones for which one of the measurement 
succeeded and the other did not. For example, if Alice's 
measurement succeeded, then she can tell Charlie the re- 
sult of her measurement, and Charlie can see whether it 
corresponds to the state that he sent. If it does not, then 
they know that an eavesdropper was present. 

A second example concerns operator discrimination 
[22-[24|]- Alice starts with the two-qubit state, \^i n ) = 
|0)|0). She sends the state through one of two black 
boxes, each black box performing an operation on the 
input state. The first black box performs an unknown, 
arbitrary single qubit rotation on the second qubit. The 
second first performs an unknown, arbitrary single qubit 
rotation on the second quibit and a Hadamard operation 
on the first, and this is followed by sending both quibits 
through a C-NOT gate, with the first qubit as the con- 
trol and the second as the target. Alice then sends the 
resulting output state to Bob, who must decide which 
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black box Alice used. Note that what is being done here 
is the discrimination between two sets of operators; the 
first black box performs an abitrary operator from the 
first sent and the second black box performs an arbitrary 
operator from the second set. 

If the input state is sent through the first black box, 
the output state that is sent to Bob is 

|¥ lout )=a|0)|0)+/3|0>|l), (3.4) 

where a and (3 are unknown. If the input state was sent 
through the second black box, Bob receives the state 

\*2out) = -^[«(|0>|0> + |1>|1» 

+/3(|0)|1) + |1)|0)). (3.5) 

The state \^>i ou t} lies in the subspace spanned by the 
vectors {|0)|0), |0)|1)} and the state \^2out) lies in the 
space spanned by {(|0)|0) + |1)|1)), (|0)|1) + |1)|0))}. That 
means that distinguishing \^i ou t) and \^2out) reduces to 
the problem of distinguishing these two supspaces. Mak- 
ing the correspondence with our previous example 

|0>|0>-|0) |1>|1>-|2) 
|0>|1>-|1) |1>|0>-|3>, (3.6) 

we see that the problem reduces to the one we have al- 
ready solved. The subspace in which \^i ou t) lies cor- 
responds to Si and the one in which \^2out) lies corre- 
sponds to 5*2. Therefore, the POVM we have already 
found will optimally distinguish (assuming that the in- 
put state is |0)|0)) through which black box the input 
state was sent. 



IV. CONCLUSION 

We have presented a POVM that optimally and unam- 
biguously discriminates between two subspaces. The con- 



struction of this POVM made use of the Jordan bases of 
the two subspaces. The results are, in fact, more general 
than what is stated in the title. They represent the com- 
plete solution to the problem of optimal unambiguous 
discrimination between mixed states of a special class, 
viz. between those states for which the spectral form 
coincides with the Jordan representation. 

We presented two applications of the measurement 
procedure, discriminating between two-particle states if 
one has only one of the particles and deciding to which 
of two sets an unknown quantum operation belongs. 

This procedure can be used to distinguish arbitrary 
mixed states by discriminating between their supports, 
but the results will not, in general, be optimal. In order 
to optimally discriminate between mixed states the struc- 
ture of the states within their supports must be taken into 
account. Based, however, on the results of this paper we 
believe that this is an extremely difficult task for density 
matrices of Rank > 2. The case of optimally discriminat- 
ing between arbitrary Rank 2 density matrices appears 
more tractable, however. How it can be accomplished is 
a problem that still remains open. 
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